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Abstract 

The renormalization group method is applied to the study of discrete dynamical 
systems. As a particular example, the Henon map is considered as applied to 
describe the transverse betatron oscillations in a cyclic accelerator or storage 
ring possessing a FODO-cell structure with a single thin sextupole. A powerful 
renormalization group method is developed that is valid correct to fourth order in 
the perturbation amplitude, and a technique for resolving the resonance structure 
of the Henon map is also presented. This calculation represents the first successful 
application of a renormalization group method to the study of discrete dynamical 
system in a unified manner capable of reducing the dynamics of the system both 
far from and close to resonances, thus preserving the symplectic symmetry of the 
original map. 
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1 Introduction and Basic Equations 



It is increasingly important to develop an improved theoretical understanding of the non- 
linear dynamics of charged particle beams in high energy accelerators and storage rings 
[H El El in El Ej • An individual particle propagating in an accelerator experiences a growth 
of amplitude of betatron oscillations in a plane transverse to the particle orbit whenever a 
perturbing force acts on it. This force may be of various origins, for example, high-order mul- 
tipole magnetic field errors, space-charge forces, beam-beam interaction force, power supply 
ripples, or other external and collective forces. 

There is a growing number of robust analytical methods used to study the effects of 
the nonlinear behavior of beams in accelerators and storage rings, ranging from classical 
perturbation theory [71IH] to the Lie algebraic approach [31101 HI]- The recently developed 
renormalization group (RG) method has been successfully applied to both continuous dy- 
namical systems fT^ IT^ ITij and maps ^3 US] that are of general interest in the physics of 
accelerators and beams. The advantage of the RG method is associated with the fact that 
it is equally powerful to study finite-dimensional, as well as continuous systems. Therefore, 
it is also useful when applied to analyze the properties of chaotic dynamical systems in both 
the stability region and the globally stochastic region in phase space [T71 ITHl ITU] . 

While the RG method is well established in applications to continuous dynamical systems, 
the present paper demonstrates that the RG method can also be applied successfully to study 
discrete dynamical systems. As a particular example, we consider the Henon map [201 1^ 
as applied to describe the transverse betatron oscillations in a cyclic accelerator or storage 
ring possessing a FODO-cell structure with a single thin sextupole. The basic equations 
and Henon transfer map used in the present analysis are summarized later in section 1, and 
in section 2 a powerful RG technique is developed that is valid correct to fourth order in 
the perturbation amplitude. A technique for resolving the resonance structure of the Henon 
map is discussed in section 3, and in section 4 illustrative numerical results are presented. 

The present analysis assumes that a certain multipole nonlinearity is concentrated at 
a single point azimuthally located at Oq. Then the one-turn transfer map describing the 
transverse betatron motion in a cyclic accelerator or storage ring can be written in the form 
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where z is the state vector 



z = 
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TZt is an orthogonal matrix 
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and F(z; 6*0) is defined by 
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Here, (X, Z) is the transverse displacement, {Px-, Pz) is the transverse canonical momentum, 
/ is the length of the multipole element, R is the mean machine radius and the 
horizontal and the vertical betatron tunes, respectively. Furthermore, V{X, Z] 6) is the 
nonlinear potential, which can be expressed in the polynomial form 



v{x,z-e) = Y: E bi'im'z^ 



rym 



(1.5) 
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where b\!^{9) are coefficients (generally functions of the azimuthal angle 9) representing the 
strength of the nonlinearity. 



2 The Henon Map 

The simplest nontrivial example of a polynomial transfer map is the so-called Henon map 
1^ . It can describe the horizontal betatron oscillations in an accelerator possessing a 
FODO-cell structure with a single thin sextupole. The two-dimensional Henon map can be 
obtained from equation in the case when the potential V{X, Z; 9) [see equation ()1.5p ] 
contains a single localized cubic nonlinearity. In explicit form it can be written as 



Xn+i = X^cost^ + i^Pn - SX^j sinu, (2.1) 
Pn+i = -Xn sinco + {Pn ~ SX^^ coscu, (2.2) 

where 

2R^ ^ ' 

In section 1, we briefly discussed the essence of the RG method and emphasized its power 
to handle a number of problems arising in the theory of continuous dynamical systems. 
In this section, we will demonstrate that the RG method can be applied successfully to 
study discrete dynamical systems, and as a particular example, we consider the Henon map. 
The latter can be further simplified by eliminating the canonical momentum variable P. 
Multiplying equation ()2.H) by coscu, multiplying equation ()2.2|) by — sincj, and summing the 
two equations, we obtain 

Xn+i cos a; — Pn+i sin a; = X„. (2.4) 
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Substitution of the recursion relation (j2.4j) into equation (|2.ip yields the second-order differ- 
ence equation for X 

CXn = Xn+i — 2Xn cosciJ + = ~eSX^ sin a;. (2.5) 

Here, e is a formal small parameter introduced for convenience to take into account the 
perturbative character of the sextupole nonlinearity. 

Next we consider an asymptotic solution of the map ()2.5|) for small e by means of the 
RG method. The naive perturbation expansion 

X. = X(°) + eX(^) + e^X^^) + . . . , (2.6) 

when substituted into equation ()2.5p . yields the perturbation equations order by order 

£Xf = 0, (2.7) 

= _5X(°)2sincj, (2.8) 
£X(2) = -25XWx«sin^, (2.9) 
£Xf = -5(X(^)^ + 2X(°)X(^)) sine, (2.10) 

£X(^) = -25(xf X(^) + X(^)xf ) sine. (2.11) 
Solving equation ()2.7|) for the zeroth-order contribution, we obtain the simple result 

X(°) = Ae'"" + c.c, Pi^^ = lAe'^'' + c.c, (2.12) 

where A is a complex constant amplitude. By virtue of equation ()2.12|) the first-order 
perturbation equation ()2.8p becomes 

£X(i) = -^(A^e^^-" + 1\Af + c.c.) sine. (2.13) 

The solution to equation ()2.13|) can be expressed as 

X(i) = -S\Af cot I + ^^lA^e^^-" + c.c, (2.14) 

where 

S, = ^— -. (2.15) 

cos 00 — COS 2uj 

To avoid resonant secular terms, we assume in addition that 

271 4:71 

u^2k7i, e 7^ h2A;7r, e 7^ h 2fc7r, (2.16) 

3 3 

where k is an integer. The properties of the Henon map in the case where the betatron tune 
is near a resonance (in particular the 1/3 resonance) will be considered in the section 3. 
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The second-order perturbation equation (|2.9p becomes 

£X^2) _ -2J^^iriu\A\^Ae'^'' - sixi^ uSSiA^e^'^"" + c.c, (2.17) 

where 

J^ = 5(^y sincu-^cot^j. (2.18) 
The solution to equation (j2.17|) can be readily expressed in the form 

) = mJ^|A|'Ae'"" + ^l^^i^sA^e^^*"" + c.c, (2.19) 

where 

^2 = — , u^{2k + l)n, uj^{2k + l)^. (2.20) 

cos OJ — COS M 2 

Continuing further with the third-order calculation, we note that 

£X^3) _ _5(^So|A|^ + Ssl^l^A^e^*"" + E4A^e^*^" + c.c.) sincj. (2.21) 
Here, the notation 

So = 5^ cot^ - + ^ sin^ u;, (2.22) 

S2 = Si (^2 sin — 5 cot — ^ sin oj + 2mjF, (2.23) 
c 

S4 = ^(5i + 4^2) sin^ (2.24) 

has been introduced. The solution to the third-order perturbation equation ()2.21|1 is found 
in a straightforward manner to be 

X'^^ = --So|A|'cot| + {B + mSiJ^smuj)\AfA^e^''^'' 

+ ^53X4^^6^'"" + C.C., (2.25) 

where 

;B = — ( ^2 sin uj — S cot — ) sin^ cu sin sin 2uj, (2.26) 

2 V 2 / o 

5 27r 47r 

53 = — , ^ — + 2A;7r, ^ — + 2A;7r. (2.27) 

cos — COS 4uj 5 

Finally, we retain terms proportional to the fundamental harmonic e**^" in the fourth-order 
perturbation equation, i.e., 

CXi^^ = {C- 2mJ^2 gij^^j |A|Me'"" + c.c, (2.28) 

where 

B+^j^sin^uj ^cot-jsinc^. (2.29) 



We obtain the secular fourth-order contribution to the fundamental harmonic in the form 



(2.30) 



where the coefficient T> is defined by 

V 



C + J-'^ cos UJ 
2 sincj 



f2.3r 



Close inspection of the naive perturbation solution starting from the second-order re- 
sult ()2.19|) shows that secular terms (proportional to n, n^, etc.) are present. To remove 
these terms, we define the renormalization transformation A A{n) by collecting all terms 
proportional to the fundamental harmonic e^^^. This gives 



A{n) 



l + it^nJ^\A\^ + t^iinV 



A. 



^2.32) 



Solving perturbatively equation ()2.32|) for A in terms of A{n), we obtain 



A 



1 — ie n!F 



A{n) +0(t 



A{n). 



^2.33) 



A discrete version of the RG equation can be defined by considering the difference 



A{n + 1) - A{n) = ie^T\AfA + eHiV 



2 



nTAm^A. 



(2.34) 



Substituting the expression for A in terms of Ain) [see equation ()2.33|) ] into the right-hand- 
side of equation ()2.34|) . we can eliminate the secular terms up to O(e^). This gives 



A(n + 1) 



1 + ie'T 



~ 2 T'^ 
A[n) +e^iV- — 



A{n] 



A{n). 



[2.35) 



This naive RG map does not preserve the symplectic symmetry of the original system and 
does not have a constant of the motion. To recover the symplectic symmetry, we regularize 
the naive RG map ()2.35p by noting that the coefficient in the square brackets multiplying 
A{n) can be exponentiated as 



A{n + 1) = A{n) exp iuj(^ A{n) ) 



where 



A{n] 



A{n] 



A{n] 



^2.36) 



^2.37) 



Although the renormalization procedure described above may seem somewhat artificial, it 
holds in all orders. By extracting a symplectic implicit map in terms of the real part and the 
argument (phase) of the renormalized amplitude v4(n), a partial proof (up to fourth order) 
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of this assertion will be presented in the next paragraph. It is clear now that the regularized 
RG map (|2.36|) possesses the obvious integral of motion 



A{n + 1) = A{n) 



(2.38) 



It is important to note that the secular terms encountered in the higher harmonics 
(e^*"^", e^*"^", etc.) can be summed to give the renormalized amplitudes, which when ex- 
pressed in terms of A{n) do not contain secular terms. This means that once the amplitude 
of the fundamental harmonic is renormalized, any problems associated with divergences in 
higher harmonics are being remedied automatically. To demonstrate this, we express the 
amplitude of the second harmonic as 



A, 



— sincu + e^{B + inSiJ^ sin u)\A\ 



A' 



which by virtue of equation (|2.33|) acquires the form 



A, 



5i 



sin + e B 



A{n] 



A\n). 



By analogy, the amplitude of the third harmonic 



A, 



3m 



sin^ uj + e^[B3 + — ^i^aJ^sin^ uj]\A\ 



A' 



can also be renormalized. The result is 



A, 



S1S2 



sin^ LO + e^Bs A{n) 



A^n). 



(2.39) 



(2.40) 



[2A1) 



[2A2) 



Proceeding in a manner similar to above, we can represent the canonical conjugate mo- 
mentum Pn according to 



Pn = iB{n)e^^^ + c.c. + higher harmonics. 



^2.43) 



where 



Bin + 1) = Bin) exp iuj(^ A(n) ) 



(2.44) 

Using now the relation ()2.4p between the canonical conjugate variables (X„,P„), we can 
express the renormalized amplitude i?(n) in terms of A{n) as 



Bin) 



iA{n) 



smcu 



e~H"+"^ -coscu 



(2.45) 



In addition, the sextupole nonlinearity shifts the closed orbit by the constant value (in 
normalized coordinates) 



tSJ 



cot ■ 



Pr. 



-X^otan ■ 



:2.46) 
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which is a common property for all odd multipole nonlinearities. 

Neglecting higher harmonics and iterating the regularized RG maps (|2.36|) and (j2.44j) . 
the renormalized solution of the Henon map can be expressed as 



X« =X,o + V27cosV^(:r;n), 
Pi'^ = Pco + ^/2JWh\J) cos^{J] n) - sin^CT; n) 



where 



(2.47) 
(2.48) 

(2.49) 
(2.50) 

smcu .... smcu 

It is evident that the integral of motion J has the form of a generalized Courant- Snyder 
invariant [231, which can be expressed as 

2 



cos u - cos [u + u{J)] 



sin [uj + uj{J)\ 



p(i) _ p 



a 



(1), 

H 



:2.51] 



It is important to emphasize that equation ()2.51|) represents a transcendental equation 
for the invariant as a function of the canonical variables {X,P), because the coefficients 
an and Ph depend on J'. Note also that the sextupole nonlinearity gives rise to a nonlinear 
tune shift uj, leading to the distortion of the invariant curves [circles in normalized phase 
space {X, P)]^ even in an approximation where only the contribution of the first harmonic 
is taken into account. Further distortions of the phase-space trajectories are introduced by 
higher harmonics. 

Taking into account all harmonics up to the fifth harmonic, we can express the renor- 
malized fourth-order solution of the Henon map in the form 



X„ = X(i) + ^ Xm cos M^iJ; n) 



M=2 



Pn = Pi^^ +Y.^M a^H^ COS Mij{ J; n) - /J^f ^ sin M^PiJ; n) 



(M) 



M=2 



The amplitudes Xm of the various harmonics are given by the expressions 

-^(Si sin a; + e^BJ'^, X4 = -e^(S3S4i7^ sina;, 



X,= '^^^{S,S2sm'u; + e'BsJ 



4 
^"5 



1 



2V2 



(2.52) 
(2.53) 

(2.54) 
(2.55) 



Furthermore, similar to equation ()2.5()|1 . the generalized a^^"* and functions can be 
expressed as 



"if' (J) 



cos UJ — COS M[lj + Uj{J')] 



smu; 



sin M[uj + uj{J)] 



smu; 



:2.56) 



tNote that an = and /3h = 1 for = 0. 
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3 Resonance Structure of the Henon Map 

The solution to the first-order perturbation equation ()2.13|) was obtained in the form ()2.14|) 
assuming that the unperturbed betatron tune v is far from the third-order resonance 3z/ = 1. 
It is important to study the properties of the Henon map near a nonhnear resonance by means 
of the RG method. In what follows, we demonstrate how the RG reduction of the Henon 
map works near the one-third resonance. A similar procedure can be performed near all 
other resonances. 

Let us expand oj according to 

2 e- ^TT 

o; = Wo + edi + e + • • • , ojq = — . (3.1 

Equation ()2.5j) can then be expressed in alternate form 

LoXn = Xn+i — 2X„ coscuo + = 2X„(costi; — coscjo) — ^SX'^ sina;. (3.2 

The perturbation expansion ()2.6|) . when substituted into equation ()3.2|) . yields the pertur 
bation equations 

£oXf = 0, (3.3 

>CoXi^) = -^(25iXf + 5X^2), (3.4 

C,X^^ = -V36,X^:^ + (I - VsS.'j + S (|x(°)^ - V3X(°)X(1)) (3,5 
Noting that 2uJo = 27t — ljq, equations ()3.3|) - ()3.5|) can be solved, yielding the result 

X(°) = Ae^^"*^ + c.c, ) = iAe^^o" + ex., (3.6 

X« = -^\Af + inG{A, A*)e'""" + c.c, (3.7 

X^2) =j^^inB+ (inC + n^pje^'^o" + c.c, (3.8 

where 

g{A,A*)=5^A+^A*'^, (3.9 

A='^S\A\\ B = ^{g*A-gA*), (3.10 

v = ]^{-5,g + sg*A*), c = -^^, (3.11 

S= (^l-y/l3^A+^-^A*-' + S^\A\^A. (3.12 

Proceeding as before, we define the renormalized amplitude by 

A{n) = A + ieng + t^(inC + n'^V). (3.13 
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Taking into account the expression 



(3.14) 



which relates the amphtude A to the renormahzed amphtude A{n), we obtain the renormal- 
ized resonant map 



where 



A{n + 1) = A{n) + ieg{n) + \iC{n) + V{n) 
X{n) =A\A{n),A*{n) 



(3.15) 
(3.16) 



and A represents C, T> or Q. 

It is important to note that the resonant shift in the closed orbit is automatically renor- 
mahzed, once the renormalization transformation A —* A{n) has been performed. The result 
is 



XM = — 



■V3 + 2e5i + 0(e^) A{ 



(3.17) 



Note that the closed orbit can be corrected up to third order (in the sextupole strength S) 
by choosing the first-order resonance detuning 5i to be 5i = \^/2. In terms of betatron 
tune, this implies 

^ (3.18) 



Az/ 



Ait 



Since the naive RG map (j3.15p does not preserve the symplectic structure of the original 
Henon map, an important step at this point consists of constructing a symplectic map 
in appropriate variables equivalent to ()3.15p . Unfortunately, the regularization procedure 
described in the previous paragraph cannot be applied to the map ()3.15p . The reason is 
that in the resonant case A{n) is no longer an integral of motion. An alternative way to 
overcome this difficulty is to represent A{n) as 

Ain) = v/Xe^^", (3.19) 
and attempt to find an (implicit) map in terms of the new variables (</?„, J'n) of the form 

^n+l = V^n + g{Jn+l, <^n) , Jn+1 = Jn + f{Jn+l, <^n)- (3.20) 

Expanding the unknown functions / and g in a perturbation series 

oo oo 

f = T.^'h, 9 = E^'9k, (3.21) 



k=l 



and substituting equation ()3.19|1 into equation ()3.15jl . we can determine / and g up to second 
order. We obtain 



gi{ip,J) = Si + —\fj cos3v?. 



(3.22) 
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/2(^, J) = 36iSJ^/^ cos 3^ - j3/2 sin 3^ + cos 6^, (3.23) 

g2{^, J) = 02 ^ ^^"^ 

- VJcos3(^- — Jsin6(^, (3.24) 
The map ()3.20|1 is symplectic provided the condition 

1^ + ^ = (3.25) 

holds. To verify equation ()3.25|) . we evaluate the determinant of its Jacobian 



det {Jm) = det 



f \ -C dj^n+l 1 _|_ f dJ^n+l 



(3.26) 



where subscripts denote differentiation with respect to the variables indicated. Taking into 
account 



diPn 
dJn 



dJn+l 










dJn+l 1 


dJn 


^ dJn 



(3.27) 



we obtain 



det (Jm) = • (3.29) 

The requirement that det (^J^m) = 1 leads to the condition ()3.25p . It is straightforward to 
verify that / and g as given by equations ()3.22j) - ()3.24|1 satisfy equation ()3.25p . 

The representation ()3.19j) of the renormalized amplitude A{n) together with ()3.2()|1 can 
be used as an alternate way to obtain the exponential form ()2.36|) of the RG map ()2.35|) . 
The expansions 



oo 

/ = Ee'V2., 9 = Y.^''92k, (3.30) 

k=l k=l 

when substituted into equation ()2.35|) . after some straightforward algebra lead to the result 

f2{v,J)=0, U{^,J) = 0, (3.31) 

92{V, J) = :FJ, g^iy^, J) = VJ\ (3.32) 

Thus, we obtain the symplectic implicit map 

V'n+l = V'n + Jn+l + C^^^i/n+l, <Jn+l = Jni (3.33) 

which is the RG map ()2.3(j|l written for the real part and the argument of the amplitude 
A{n). 
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4 Numerical Results 



In this section we present illustrative numerical results revealing the stability properties of 
the Henon map. We use the analytical renormalized solution expressed by equations ()2.52|) 
- fl2.5(i|l to construct the phase portrait of the Henon map near the third-order resonance 
with V = 0.323 (also u = 0.34525), near the fourth-order resonance with u = 0.24, and near 
the fifth-order resonance with u = 0.19. All calculations are performed for a relatively large 
value of the sextupole strength corresponding to S = 0.1. 

In Figure 1 the phase portrait of the Henon map near the third-order resonance with 
u = 0.323 is depicted. It shows the stability region and the invariant curves for different 
values of the modulus squared of the renormalized amplitude A [see equation (|2.38j) ]. As the 
value of the invariant J' increases, it reaches a value J^max above which the phase trajectories 
begin to intersect. This is due to the fact that the perturbation renormalization technique, 
valid for unperturbed betatron tunes sufficiently far from resonances, does not work well 
in this region, and the reduction procedure near resonances developed in section 3 should 
be employed. However, the quantity Jimax is closely related to the dynamic aperture. For 
the case where z/ = 0.323 it is approximately JTmax = 3.85. Figure 2 represents the phase 
portrait of the Henon map for u = 0.34525. The value of J^max is the same as in the case 
where u = 0.323. 

Figures 3 and 4 show the phase portrait of the Henon map near the fourth-order and the 
fifth-order resonances where u = 0.24 and u = 0.19, respectively. The dynamic aperture in 
the case of the fourth-order resonance is approximately J^max = 19.01, while in the case of 
the fifth-order resonance it is J'max = 28.5. 

5 Conclusions 

While the renormalization group method is well established in applications to continuous 
dynamical systems, the present paper demonstrates that the renormalization group method 
can also be applied successfully to study discrete dynamical systems. As a particular exam- 
ple, we considered the Henon map as applied to describe the transverse betatron oscillations 
in a cyclic accelerator or storage ring possessing a FODO-cell structure with a single thin 
sextupole. The basic equations and Henon transfer map used in the present analysis were 
summarized in section 1, and in section 2 a powerful renormalization group technique was 
developed that is valid correct to fourth order in the perturbation amplitude. A technique 
for resolving the resonance structure of the Henon map was discussed in section 3, and in 
section 4 illustrative numerical results were presented. To the best of our knowledge, the 
present calculation represents the first successful application of a powerful renormalization 
group method to the study of discrete dynamical systems in a unified manner. In section 3 
it was shown that the renormalization group map can be further expressed in terms of an 
implicit symplectic map in both cases, far from and close to resonances. Further applications 
to discrete dynamical systems will include generic polynomial transfer maps, the standard 
Chirikov- Taylor map, etc. 
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Figure 1: Phase portrait of the Henon map obtained from equations ()2.52j) - ()2.5(j|l near the 
third-order resonance with v = 0.323. Here, takes values ranging from J' = 0.15 (inner 
contour) to J' = Jmax = 3.85 (outer contour). 
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Figure 2: Phase portrait of the Henon map obtained from equations ()2.52j) - ()2.5(j|l near the 
third-order resonance with v = 0.34525. Here, takes values ranging from = 0.15 (inner 
contour) to JT" = J^max = 3.85 (outer contour). 
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Figure 3: Phase portrait of the Henon map obtained from equations ()2.52j) - ()2.5(j|l near the 
fourth-order resonance with v = 0.24. Here, takes values ranging from J' = 0.5 (inner 
contour) to JT" = J^max = 19.01 (outer contour). 
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Figure 4: Phase portrait of the Henon map obtained from equations ()2.52j) ■ 
the fifth-order resonance with v = 0.19. Here, takes values ranging from J' 
contour) to JT" = J^max = 28.5 (outer contour). 
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